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Abstract. A method is developed for validating model-based estimates of atmospheric 
moisture and ground temperature using satellite data. The approach relates errors in 
estimates of clear-sky longwave fluxes at the top of the Earth-atmosphere system lo errors 
in geophysical parameters. The fluxes include clear-sky outgoing longwave radiation 
(CLR) and radiative flux in the window region between 8 and 12 /x:n (RndWn). The 
approach capitalizes on the availability of satellite estimates of CLR and RadWn and 
other auxiliary satellite data, and multiple global four-dimensional data assimilation 
(4-DDA) products. The basic methodology employs off-line forward radiative transfer 
calculations to generate synthetic clear-sky longwave fluxes from two different 4-DDA 
data sets. Simple linear regression is used to relate the clear-sky longwave flux 
discrepancies to discrepancies in ground temperature (57^) and broad-layer integrated 
atmospheric prccipitablc water (5/;w). The slopes of Lite regression lines delinc sensitivity 
parameters which can be exploited to help interpret mismatches between satellite 
observations and model-based estimates of clear-sky longwave fluxes, bor illustration we 
analyz.c the discrepancies in the clear-sky longwave fluxes between an early 
implementation of the Goddard Earth Observing System Data Assimilation System 
(GEOS2) and a recent operational version of the European Centre for Medium-Range 
Weather Forecasts data assimilation system. The analysis of the synthetic clear-sky flux 
data shows that simple linear regression employing 57’,, and broad layer 5 /mv provides a 
good approximation to the full radiative transfer calculations, typically explaining more 
than 90% of the 6 hourly variance in the flux differences. These simple regression 
relations can be inverted to “retrieve” the errors in the geophysical parameters. 
Uncertainties (normalized by standard deviation) in the monthly mean retrieved 
parameters range from 7% for 8T g to ~~20% for the lower tropospheric moisture between 
500 hPa and surface. The regression relationships developed from the synthetic flux data, 
together with CLR and RadWn observed with the Clouds and Earth Radiant Energy 
System instrument, arc used to assess the quality of the GEOS2 I g and/nv. Results 
showed that the GEOS2 T g is too cold over land, and pw in upper layers is loo high over 
the tropical oceans and too low in the lower atmosphere. 


1. Introduction 

While much progress has been made to improve the climate 
characteristics of general circulation models (GCMs), the hy- 
drological cycle stands out as a major component of the Earth- 
Atmosphere system which is still poorly modeled |c.g., Gates ct 
ol 1999; Latt ct <*/., 1995]. Major advances in modeling the 
hydrological cycle arc hampered by inadequate and/or incom- 
plete measurements of such quantities as precipitation, latent 
heating, clouds, atmospheric and soil moisture, and ground 
temperature. Data obtained by methods of four-dimensional 
data assimilation (4-DDA) suffer from errors in these quanti- 
ties, both as a result of errors in the assimilating models and a 
lack of observations lo directly constrain the hydrological cycle. 

This paper is not subject lo U.S. copyright. Published in 2000 by the 
American Geophysical Union. 
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[n this study we focus on methods to verify model-based 
(including 4-DDA) moisture profiles (r/) and ground temper- 
ature (T t/ ) using satellite data. We focus on these two param- 
eters because they arc highly dependent on the parameteriza- 
tions of sub-grid-scale processes, and they arc difficult to 
validate because there are Tew reliable observations of these 
quantities with global coverage. The quality of the ground 
temperature directly reflects on the quality of the land surface 
formulation, while the moisture profiles arc importa it con- 
straints on the behavior of the boundary layer and convection 
schemes. The methodology described here was developed as 
the result of efforts to validate these quantities in a global data 
assimilation system, but the methodology should also prove 
useful for assessing systematic errors of global almt spheric 
models. 

The basic methodology is as follows; First, sensitivity param- 
eters arc obtained by applying off-line forward radialiv ; trans- 
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fer calculations to two different assimilated data products to 
relate differences in dear-sky longwave fluxes at the top of the 
Earth-atmosphere system U) errors in geophysical parameters. 
The clcar-sky longwave fluxes indude dear-sky outgoing long- 
wave radiation (CLR) and radiative flux in the window region 
between 8 and 12 pm (RadWn). The assimilated data sets 
should provide physically reasonable and global estimates of 
the geophysical fields. Any discrepancies in the fields from the 
two different systems measure our uncertainly and provide a 
tool for assessing the sensitivity of the CLR to discrepancies in 
these fields. 'Hie sensitivity parameters together with satellite 
estimates of clcar-sky longwave fluxes (and other auxiliary sat- 
ellite data) arc then used to relate the clcar-sky longwave 
fluxes to errors in the geophysical fields. For purpose of illus- 
tration, sensitivity parameters arc obtained by comparing 
ground temperature anti atmospheric temperature anil mois- 
ture profiles from an early implementation of version 2 of the 
Goddard Earth Observing System (GEOS2) data assimilation 
system [Data Assimilation Office (DAO), 1996] with those 
quantities from a recent operational version of the European 
Centre for Medium-Range Weather Forecasts (ECMWF) data 
assimilation system [Courtier cl ai, 1998]. 

The key satellite data used for this study are the CLR and 
RadWn estimates from the Clouds and Earth Radiant Energy 
System (CERES) [Wiclicki ct ai , 1996] instrument onboard the 
Tropical Rainfall Measurement Mission (TRMM). 

This approach to validation is viewed as a compromise be- 
tween a simple comparison of model-generated and observed 
CLR and RadWn and full retrieval algorithms that attempt to 
exploit information in multiple spectral bands to obtain esti- 
mates of the geophysical quantities. The method provides in- 
sight into the nature of the errors in the model-based CLR and 
RadWn consistent with the broadband nature of the CLR and 
RadWn satellite measurements. 

Sections 2 describes the satellite data and radiative transfer 
model used for the forward calculation. The methodology is 
described in section 3. The sensitivity analysis, based on the 
GEOS2 and ECMWF assimilated data, is described in section 
4. A comparison with the CERES/TRMM clear CLR and 
RadWn is presented in section 5. Discussion and conclusions 
are given in section 6. 

2. Satellite Data and Radiative Transfer Model 

The primary quantity analyzed in this study is the CLR and 
RadWn. The availability of high-quality satellite estimates of 
CLR and RadWn (sec below) has made these important quan- 
tities for validating global models. Relating and understanding 
differences between model-computed and observed CLR and 
RadWn is complicated by a number of factors. Some of the 
differences arise From sampling differences and fundamental 
differences between satellite measurements and model “grid- 
scale" fields. In this study we focus on the problem of relating 
the errors (computed minus observed) in CLR and RadWn to 
deficiencies in the model-based estimates of the radialivcly 
important geophysical quantities (c.g., atmospheric moisture 
and temperature and ground temperature). The methodology 
is developed with synthetic flux data to circumvent the sam- 
pling and representativeness problems. These issues arc be- 
yond the scope of this initial study, though they will need to be 
addressed when developing quantitative estimates of errors in 
the geophysical quantities based on satellite measurements. 
We note that in this study the model-based CLR and RadWn 


arc computed in an ofr-finc mode (see below). Many institu- 
tions now routinely compute CLR on-line during the course of 
model integrations. 

In the following we describe the CLR and RadWn data, 
various other auxiliary data sets used for this study, I nd the 
radiative transfer scheme used for the off-line calculi lion of 
CLR and RadWn. 

2.1. CLR and RadWn Satellite Data 

While the sensitivity parameters will he developed w ih syn- 
thetic flux data (see next section), our application of th*: meth- 
odology to validating GEOS2 (sec section 5) will rcqi ire sat- 
ellite observations as input. For the latter we have o Unined 
CERES/TRMM data [Wiclicki ct al. , 1996]. The qualil; of the 
CERES data is comparable to the quality of the Earth Radi- 
ation Budget Experiment (ERBE) data for instantaneous ra- 
diance, (luxes, ami scene lype. Generally, radiance uncertain- 
ties arc at the 1% level or less. Some differences tetween 
CERES/TRMM and ERBE-ERBS are as follows: the field of 
view resolution, the spectral response of the instrumci-ls, and 
the tropical only coverage of TRMM [Wiclicki ct ai, 1 )96]. 

2.2. Radiative Transfer Scheme 

For the sensitivity calculations we employ the rrdiative 
transfer scheme developed by Chou and Suarez [1994]. This 
scheme is used in the GEOS2 model. The longwave (LW) 
calculation has nine bands (band 10 is added to compute flux 
reduction due to N 2 0 in the 15 /xm region). The transmission 
and absorption oMI^O, 1EO continuum, C0 2 , O*. CEC, CII.,, 
and N 2 0 is modeled using & distribution. In the LW, multipli- 
cation approximation for bands arc assumed. 1 he LW scheme 
compares well with detailed line-by-line calculations. The root 
mean square (rms) errors of CLR arc between I and 3 Wm 
[Chou and Kou van's, 1991; Chou ct ai, 1995; Kratz ct al., 1998]. 
In addition, the code has been participated in the Intercom- 
parison of Radinlion Codes in Climate Models (ICECCM) 
[EUin^son ct id., I99I|. 

Surface cmissivity depends on surface property or vegetated 
surface. For a given surface type, its values are also different 
for different wavelengths in the LW spectrum. The global 
distribution of surface type and the associated cmissivity that 
we used were taken from CERES. Hie Advanced Spaccbomc 
Thermal Emission and Reflection Radiometer created ail eas- 
ily accessible data set on the basis of extensive measurements 
of the spectral reflectances of surface material in the 2 -16 pm 
region by Salisbury and D'Aiia [1992]. Wilber ct al. [1999] de- 
rived spectrum cmissivity from the spectrum reflectance on the 
basis of scene type. The scene type was determined by using a 
1 km map of the International Geosphere Biosphere Program 
scene types as supplied by the U.S. Geological Survey (f JSGS). 
A scene type of tundra was added to separate it froiT desert 
resulting in 18 surface types. We calculate the fraction jf each 
scene type in a 2° X 2.5° longitude and latitude box fi Jin the 
1/6° equal angle data. Radiative fluxes arc computed in :ach of 
the 2" x 2.5° boxes several times on the basis of the number of 
scene types in a box (maximum is 18). The final mean r; diativc 
flux is an area mean of the fractional type. Surface cmissivity 
effect on CLR is large over the desert regions, where iff-Iinc 
computations show it can lie as high as 5-8 Wm 2 mostly 
from window bands). 

2.3. Ozone 

The analyzed ozone fields from the Goddard Earth Observ- 
ing System (GEOS) ozone data assimilation system (DAS) 
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[Riishpjgaard ct al 1999; Stajner ct al., 1999] arc used in this 
study. They arc obtained by assimilating total column ozone 
observations from the Total Ozone Mapping Spectrometer 
(TOMS) and ozone profile observations from the solar baek- 
scatlcrcd ultraviolct/2 (SBUV/2) instruments into an off-line 
transport model driven by the GEOS-DAS assimilated winds. 
While the statistical analysis of ozone is performed at every 
model lime step (15 min) using a global physical space-based 
analysis scheme, the analyzed ozone is written out at the same 
frequency and resolution as the other GEOS-DAS fields. The 
high quality of the ozone fields is illustrated by two examples 
from their validation for January 1998. The average rms dif- 
ference between TOMS observations and total column ozone 
forecast is about 1 1 Dobson units or 3.7% of the average total 
column ozone. The comparison of analyzed ozone profiles with 
independent ozone observations measured by the Halogen Oe- 
cullation Experiment (HALOE) onboard the NASA Upper 
Atmosphere Research Satellite shows that the mean profiles 
agree within 3.2% between the pressures of 30 and 1 hPa. In 
this region the HALOE measurements agree within 5% with 
ozoncsondc, lidar, balloon, rockctsonde, and other satellites 
ozone measurements [Bruit! ct al 1996]. 

While the ozone contribution to CLR is relatively small, the 
accuracy of its distribution is critical for the CLR computation. 
In particular, the high ozone concentrations in the stratosphere 
and its presence in the water vapor window region are impor- 
tant for the clcar-sky flux. 


3. Methodology 

To help understand the methodology, we start with the ra- 
diative transfer equations for clcar-sky conditions in the LW 
region of the spectrum given 



f n f" </r r (u - //') 

+ vDjLnu’))--^ — itu' du t (i) 

Jn -'ll 

where v is frequency, B u is the Planck parameter, e w is surface 
cmissivity, t„ is the transmission parameter at frequency u, u is 
the path length, and u, is the total path length. In (1) the first 
terms on the right-hand side (RHS) is the contributions from 
the surface. The second term is the downwclling atmospheric 
radiation reflected from the surface. This term tends to be 
small compared to the other terms. The third term is the 
contribution from the atmosphere. 

The first step is to develop sensitivity parameters relating 
CLR differences to differences in the geophysical quantities. 
This is carried out here using data from two different four- 
dimensional data assimilation systems. The simulated radi- 
ances arc computed using the T fJ% and q and temperature (T) 
profiles as input to the radiative transfer calculations using the 
scheme described in the previous section. The aim of this step 
is to try to understand which aspect of the input data is most 
sensitive to the CLR and how the sensitivity varies over the 
globe. 


We consider a change in the CLR in terms of the following 
linearization: 
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where 5r/ ( and arc the differences in the two moisture and 
temperature products in layer i. We consider the par ial de- 
rivatives in (2) as sensitivity parameters relating change; in the 
geophysical quantities to changes in the CLR. The validity of 
this approximation will be examined in the following sci lion by 
comparing results from the full radiative transfer calculations 
(referred to as the "true” values in the following discissions) 
with results from "best lit” linear regression equations based 
on (2). 


4. Sensitivity Analysis 

In this section we shall compute the sensitivity parameters 
defined previously using output from two different l-DDA 
systems (GEOS2 and ECMWI 7 ). Note that since we arc only 
interested in differences, there is no need in these calculations 
to assume one or the other 4-DDA systems is concct. In 
section 5 we will use the sensitivity parameters together with 
satellite data to validate the GEOS2 system. 

GEOS2 represents a major upgrade to the baseline CEOS I 
system employed in the NASA first rcanalysis effort \Schnhcrt 
ct al., 1993]. 'Hie final version of CEOS 2 includes a physical- 
space three-dimensional variational analysis algorithm (Physi- 
cal-Space Statistical Analysis System (PSAS) [Coin t cl a!., 
199S]) and numerous improvements to the general circulation 
model. The model improvements include a more accurate dy- 
namics [Takacs and Suarez , 1996], a gravity wave drag scheme 
[Zhou cl al., 1996], improved diagnostic clouds, the mosaic 
land surface scheme of Kostcr and Suarez [1994], a Uvcl 2.5 
moist turbulence scheme f Hclfand ct al ., 1999], and n;w SW 
and LW radiation code [Chou and Suarez, 1994]. 'Hie system 
also includes the capability to assimilate TO VS moisture and 
Special Sensing Microwavc/Imagcr (SSM/I) total prcc pit able 
water. ^ U t j I '\ i 4 } 

The version of GEOS2 used here was run with 70 /erlical 
layers extending to 0.01 hPa and with a spatial rcsolutn n of T 
x 2.5° latitude-longitude. Tire T ff output was saved wery 3 
hours, while the upper air T and q were saved every t hours. 
Tie GEOS2 data arc compared with operational analy; cs gen- 
erated by the ECM WF for January 1998. These data ar : avail- 
able every 6 hours and have been interpolated to a 2° x 2.5° 
grid to be consistent with the GEOS2 data. 

As described in the previous section, T ff , T, and q frem both 
analyses arc input to a radiative transfer model to compute the 
CLR and RndWn. The primary interest for this study i; in the 
sensitivity to T tf and q. Our initial analysis showed that T 
differences between CE.OS2 and ECM WE arc small, and the 
effect on the computation of CLR is around ±1.5 WnT 2 . 

In the following, CLR computed by using all GEOS2 data is 
refer red to as CLRGEOS2, and the CLR computed using all 
ECMWEdata is referred to as CLRECMWF. Sensitivities are 
estimated by recomputing the CLR after replacing selected 
input fields. For example, to estimate r)CLR Ji)T, n we compute 
the CLRGEOS2 as before but replacing the GEOS2 T ff with 
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ihc T g from Ihc ECMWF data. We will refer to this as 
CLRGEOS2( 57^), indicating the calculation is carried out 
with the GEOS2 ground temperature perturbed by the amount 
5T f , = 7„(ECMWF) - 7^(GEOS2). The difference 
6CLR(S7 tf ) = CLRGEOS2(flT tf ) - CLRGEOS2 which di- 
vided by ST g is a measure of the sensitivity of CLR to ground 
temperature. For example, the change in CLR due to differ- 
ences in T 9 is then denoted by 5CLR(57^). An estimate of the 
sensitivity of CLR to T g changes is 

rtCLR SCLR(57\) 

0T, “ ST, ' 

Since the sensitivity of CLR to moisture varies with height, 
sensitivity ratios arc estimated separately for different layers in 
the atmosphere. For these calculations we used ECMWF data 
as the basis for determining the moisture sensitivity since we 
felt that the ECMWF moisture values were more reasonable 
especially at upper levels. In principle, it should not matter 
which one is used (GEOS2 or ECMWF), since we arc only 
after the sensitivity. In fact, tests showed similar results using 
GEOS2, though we felt the analysis using ECMWF provided 
results that were somewhat easier to interpret for the moisture. 
The layer-integrated moisture is referred to as the prccipitablc 
water (piv). For example, we compute SCLRfSpw,,) = 
CLRECMWF(fipiv,,) - CLRECMWF. Here 15piv, f is equal to 
pw from ECMWF minus pw from GEOS2, where pw h is the 
moisture in the layer between 200 and 500 hPn. Similar calcu- 
lations arc done for moisture between 200 and 700 liPa 
( /'"Goo), between 200 and 900 liPa (pw 9n() ), between 200 liPa 
and surface (pw T/c ), between 500 and surface (pw f ) t and the 
full moisture profile (tpw). The sensitivities to pw arc given by 
equations analogous to (3). For example, 

dCLR SCLR(SpwJ 

< lpw h RpWh ^ 

estimates Ihc sensitivity of CLR to thepw between 200 and 500 
hPa. 

We use the period January 1998 to compute the sensitivity 
parameters, since we have in-house operational analyses from 
ECMWF for this time period. Unless noted otherwise, all 
calculations arc based on 6 hourly data on the GEOS2 2° X 
2.5° grid. We estimate the sensitivities locally using simple 
regression to relate the 6 hourly differences in CLR (or 5CLR) 
to differences (8x) in the input quantities. For a single pre- 
dictor the regression equation takes the form 

5CLR(5jr) = + r. t (5) 

where x is the quantity being varied (for example, T, r pw h , or 
pw,), is the sensitivity parameter to be estimated, and e is 
the component of the 5CLR not explained by fix. 

Before computing the sensitivity parameters, we first exam- 
ine the differences in the CLR from the two assimilated data 
sets and look in some detail at the various hand contribu- 
tions to the total CLR. Figure I shows the mean 8CLR = 
CLRGEOS2 - CLRECMWF for January 1998 computed 
from the 6 hourly data. Over lane! the differences can be ns 
large as 15-20 WrrT 2 and over the oceans as large as 20-25 
Wm“ 2 . Figures lb and 1c show 5CLR(57 y ) and 5CLR(5^), 
respectively. The pattern and magnitude of the 6CLR distri- 
bution over land is very close to SCLRfST,,), while over 


Clear Sky OLR Dif (Wm"’) January 1998 
(Contour intervals (SWm* 1 , starting ot -5) are for netogive ofues) 



6CLR(6q) 



Figure I. Monthly mean 5CLR (Wm -2 ) based >n the 
GEOS2data and the- ECMWF data: (a) differences ind iced by 
<57 ’ 87, and 57; (l>) differences induced by 57^; ; nd (c) 

differences induced by 87. Contours arc for negative values; 
ihc intervals are the same as the shaded values. 


oceans, 5CLR is very similar to 5CLR(57). This indica cs that 
for these two assimilated data sets the largest discrepancies in 
CLR over the land poleward of about 40° latitude arc due to 
ground temperature differences, while the largest discrepan- 
cies in CLR over the oceans arc due to differences in the 
moisture profiles. 

To better understand the sources of the CLR from the 
atmosphere and land, we look in more detail at the various 
band contributions to the CLR (Figures 2 and 3). The gaseous 
absorption included in the radiative transfer code [Chou and 
Suarez, 1994] arc H 2 0, C0 2 , 0 3 , and trace gases (see section 
2.2). fn the following we will focus on the water vapor effect. 
Column I in Figure 2 shows the CLR computed from the 
ECMWF data in each of the nine bands. The first three bands 
arc the contributions from the water vapor rotational bands. 
The next three arc the contributions from the 800 to 1215 
cm -1 window region, and the last three represent contribu- 
tions from water vapor vibration and rotation bands. Column 2 
shows the emission from the surface (c„7„). Column 3 is the 
contribution to the CLR from the first two terms on the RI IS 
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ECWWF Dato (Wm’’) Jan 1 00Z 1998 


1 : 0 - 340 

2: 340 - 540 

3: 540 - 300 

4: 800 - 980 

5: 980 - 1100 


6:1100 - 1215 
7:1215 - 1350 
8:1380 - 1900 
9:1900 - 3000 


Fieurc 2 Contributions to CLR (Wnf 2 ) from different parts of the atmosphere and different bands. From 
S?S£ t,2Scs from bald 1 to 'band 9. the wave number for each band ^ 

t ] ic figure Value cT is surface emission, Surf (surface contribution) is the (irst two terms of the RHS of (1) 
from the atmosphere) is the third term of the R1 IS of (l). Contour intervals arc 10 

Wm" 2 except the last one in a column, which is 5 Wm 


of (1), and column 4 shows the contribution to the CLR from 
the atmosphere (third term on the RHS of (1))- 
The strength of surface emission follows the Planck function 
which peaks in the third band (540-800 cm '). The intensity 
from the other bands decreases gradually toward both sides. 
Note that the band widths arc not uniform. The atmospheric 
transmission function strongly modifies the surface emission, 
especially in bands with strong gaseous absorption. The degree 
of the atmospheric effect depends on the opacity of tire atmo- 
sphere. Comparing, in Figure 2, the surface emission with the 
amount reaching the top of the atmosphere, we find a signif- 
icant greenhouse effect due to gaseous absorption. Tlic green- 
house effect is very strong in the water vapor rotational blind 
ranging from zero to 800 cm" 1 (which include the 15 jim C0 2 
band) and the vibrational anti rotational band ranging from 


1215 to 3000 cm’ The largest contributions to CLR arc from 
bands 2 to 4, ranging from 340 to 980 cm *, and most of the 
contributions arc from the atmosphere because of the opacity 
of the water vapor rotational bands and C0 2 bands. For ex- 
ample, the surface has the largest emission in band 3, but the 
surface contribution to CLR is relatively small due to the 
opacity of the atmosphere. In particular, over the tropical 
oceans the surface contribution to the CLR from this >and is 
about 5%, while the rest is front the atmosphere (see Figure 3). 
Figure 3 includes the effective pressure level of the cc ntribu- 
tions from the atmosphere averaged over the month. This is 
determined by finding the peak of the weighting functic n com- 
puted from the mean transmission functions. In bands and 3, 
most contributions over land arc from layers close to airfacc 
(lowest 300 hPa), while most contributions over the 'ropical 
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ECMWF Data (Wrrr*) Jan 1 00Z 1998 


0 - 

340 



340 - 

540 

6:1 100 - 

1215 

540 - 

800 

7:1215 - 

1380 

800 - 

980 

8:1380 - 

1900 

980 - 

noo 

9:1900 - 

3000 


Figure 3. Percentage contributions to CLR (%) from Surf ami Atm arid mean peak of the weighting 
function (hPa). 


oceans arc from layers ranging from 200 to 500 hPa, where 
water vapor concentration is high. Away from the tropics the 
opacity decreases, and the layer of largest emission is lowered 
to between 500 and 700 hPa. In the tropics the atmosphere 
contributes more than the surface in all bands. The largest 
contributions are from between 200 to 500 hPa and some arc 
from 500 to 700 hPa depending on the opacity of tire atmo- 
sphere. 

Tn summary, the major contribution to CLR over the trop- 
ical oceans (30°S to 10°N) is from the atmosphere between 200 
and 500 hPa. In the subtropics the contribution tends to be 
from lower in the atmosphere between 500 and 700 hPa. Over 
land, for the Northern Hemisphere winter, most of the contri- 
bution to CLR is from the surface (from window bands) and 
from the lowest (300 hPa) layer of the atmosphere. 

Figure 4 shows the regression results for T tJ (land only) for 
January 1998 based on 6 hourly data. We note that sea surface 
temperature (SST) is specified from observations in both anal- 


yses, and we do not consider any differences in the SS T data. 
The regression is carried out on the CLR differences resulting 
from perturbing only the ground temperature using tl e radi- 
ative transfer equation (1). The extent to which the regression 
explains the results from the full radiative transfer calculation 
is thus a measure of the adequacy of the linear approximation 
for the e(Tec! of T, r The estimated sensitivity parameter 
i'jCLR/JT,, ** a r (Figure 4a) is highest in the middle and high 
latitudes (ranging from 1.2 to 1.6 Wm' 2 per 1°C change in T ff ) 
and lowest over the tropical land masses. This is consistent with 
our previous analysis of the various contributions to the CLR 
(Figure 2) and reflects the latitudinal dependence of the at- 
mospheric moisture content. The explained variance is shown 
in Figure 4b. Remarkably, the regression line (5) explains more 
than 90% of the variance over most areas (Figure 4b). These 
results suggest that the sensitivity of CLR to changes in T ff can 
be reasonably estimated by the simple linear approximation (5) 
of the radiative transfer equation. 
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ECMWF and GE0S2 January 1998 
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Figure 4. Estimation of ST, (rased on linear fit: (a) uT ;l (Wm 2 °C '), contour intervals arc 2 Wm Z o C 
(b) total variance explained; (c) monthly mean f>T ff ( D C); (d) estimated 5CLR(fi/ ,,) (Wm 2 ) based on Figures 
4a and 4c; and (c) errors (°C) in estimation of 5 T g based on Figure 4d and (f) 95% fiducial interval. Contours 
arc for negative values in Figures 4c~4f; the intervals arc the same as the shaded values. 


We arc interested in the ability of the regression to represent 
the systematic difference in the CLR fields. Figure 4c shows 
the mean ground temperature difference between the two 
data sets for January 1998. The mean temperature differences 
arc quite large over the cold continents accounting for most of 
the differences in CLR over those regions. Figure 4d shows the 
mean CLR difference from the regression equation evaluated 
for the lime mean bT fJ . Comparison with Figure lb shows that 
the linear approximation (the regression line) with just ST t/ 
provides a good approximation It) the time mean, differences in 
CLR over land. The reader is reminded that the results in 
Figure lb arc computed with the radiative transfer equation 
and include the effects of <7, 7\ and T, r We shall see in the next 


section that the differences in T {J for the most part reflect a 
cold bias in this version of GEOS2. 

We next obtain an estimate of 5 T g by inverting the regres- 
sion equation (5): i.c., by solving for 5 T ff using the estimate of 
a r . For brcvily we refer to this estimate of *)T tJ as the “re- 
trieved'' value. Figure 4e shows the difference between our 
retrieved 5 T {t and the actual values shown in Figure 4c. The 
magnitude of the errors in ST g arc in most regions less than 
PC. Figure 4f shows that the 95% fiducial intervals tor the 
inverse regression {Draper and Smith , 1 981] or retrieved values 
are typically Note that here we only show one side of the 

two-sided confidence interval corresponding to the sign of the 
actual error in Figure 4c. 


1 
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Figure 5. Same as in Figure 4 except for units are (Win ~ 2 min' l ). Contour intervals arc for values 

opposite to the shaded values. 


F5 Figure 5 is the same as Figure 4 except for the mean mois- 

ture difference in the layer between 200 and 500 hPa (5piv /( ). 
Figure 5a shows the sensitivity parameter (a riv J determined 
from the regression 

5CLR(fy;w, t ) = a imh Spw k + e. (6) 

The spatial pattern of the sensitivity to upper level moisture is 
largely zonally symmetric, with the lowest sensitivity (<3 
Wm -2 /mm prccipitablc water in the 200-500 hPa layer) in the 
convective regions of the tropics, and increasing sensitivity 
away from the tropics. A region of enhanced sensitivity (>9 
Wm~ 2 /mm prccipitablc water in the 200-500 hPa layer) occurs 
over the subtropical regions of the North Pacific. We note that 
bands 2, 7, and 8 (Figure 2) show a tendency for enhanced 
contributions to the CLR from the subtropics especially over 
the North Pacific, suggesting the increased sensitivity is coming 
largely from these bands. Figure 5b shows that the regression 
line explains more than 90% of the variance in the CLR over 


most of the globe. Scattered regions of the tropics and sub- 
tropics, the cxtratropical storm tracks, and the Himalayas, ex- 
plained variance less than 80%, With those exceptions the 
linear approximation to the full radiative transfer calculations 
for the impact of the upper level moisture appears to le quite 
good. Figure 5c shows the mean difference in the upp :r level 
moisture (itynv,,). The differences tend to be n :gativo 
(GJ20S2 is wetter) and largest just outside the tropi :s. The 
results of the regression for the time mean 5CLR(5/;u', ) (Fig- 
ure 5d) again show a good approximation to the full call ulation 
(Figure le), suggesting that much of the systematic dif crcncc 
between GHOS2 and ECMWF CLR over the subtropics arc a 
result or differences in the upper level moisture. The results of 
the inverse calculation (Figure 5e) show that the magnitude of 
(he errors in the retrieval of is in most regions less than 
0.5 mm. Figure 5T displays the 95% fiducial intervals on the 
retrieved values. The typical values range from 0.3 to C.5 mm. 

Figure 6 is the same as Figure 5 except for the mean mois- F6 







I 


TGI 

Docket #-2118068 


WU ET AL: ASSESSING GROUND TEMPERATURE AND ATMOSPHERIC MOISTURE 
ECMWF and GE0S2 January 1998 


a 


pw L 


90 


9Q 

60 

1 

1 

‘.3 

Ti. 

60 

30 


30 

0- 


0 

-30- 


-30 



-60 

i f 1 1 j 

-60 

on * 


c 

60 120 180 240 300 3t 

-90 

>0 ( 


*mmmm 

— 1.5 -0.5 


-4 


-2 


-8 


Var Expi 



60 120 180 240 300 360 

mmmmm 

0.8 0.9 





2 4 6 8 

95% fiducial interval 



-8 


-2 


Figure (%. Same ns in Figure 4 except Tor fi/Jiv,: units are (Win" 2 mm" 1 ). Contour intervals are for values 
opposite to the sliadcti values. 


turc difference in the layer between 500 hPa and the surface 
(5/nv,). Figure 6a shows that the sensitivity parameter (a /Mi>/ ) 
is about an order of magnitude smaller than the upper level 
sensitivity parameter (cr /f ^). The greatest sensitivity (magni- 
tude > I to 2 Wm -2 /mm prccipitablc water) occurs at high 
latitudes in regions that are probably ice covered. Figure 6b 
shows the variance explains by the regression line which ex- 
plains more than 90% of the variance in the CLR over most of 
the subtropics. Generally, less than 80% of the variance is 
explained over the cold continents and the polar regions. Fig- 
ure 6c shows the mean difference in the lower-level moisture 
(fynvj). TIic differences tend to be positive (GEOS2 is drier) 
throughout much of the tropics and subtropics, with the largest 
values occurring off the west coasts of Africa and South Amer- 
ica in regions characterized by the presence of low-level stratus 
clouds. Tli ere arc some regions with negative values in (lie 
Northern Hemisphere subtropics that coincide with large neg- 
ative values at upper levels (compare Figure 5c). The results of 


the regression for the time mean 5CLR((5/;iv / ) is shown in 
Figure 6d. 'Fliis shows that some of the negative values along 
the west coasts of Africa, South America, and Africa in the full 
calculation (Figures la and lc) arc due to the differences in the 
low-level moisture. Figure 6c shows that the magnitude of the 
errors in the retrieval of ftpw t is in most regions less than 2 
mm. 'Ilie 95% fiducial intervals (Figure 6f) are in the range 
±2-4 mm with the largest uncertainties in the Southern Hemi- 
sphere and the north polar region. 

We next look more closely at the vertical dependence of the 
sensitivity of CER to moisture. The panels on the left-hand 
side of Figure 7 show the differences between the ECMWF F 7 
and die GEOS2 time mean moisture. The sensitivity parame- 
ters arc the slopes of the regression lines through die scatter- 
plots of 5CLR versus 8/or (panels on die right-hand side of 
Figure 7). These arc computed, in this ease, from the monthly 
mean fields, and each point in the scattcrplots represents a 
different grid point (ocean only). This clearly shows that 5CLR 
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is most sensitive to 8pw in the top layers, flic sensitivity de- 
creases the farther the layer is from the lop of the atmosphere. 
Thus the moisture in the layer dose to the surface tends to 
have little effect on CLR, consistent with our previous analysis 
of the contributions to the total CLR (Figure 2). Also, above 
200 hPa the amount of moisture is too small to have an np- 
rrwiahfe impact on the CLR despite the strong sensitivity. 
This is illustrated in Figure 8, which shows the relative contri- 
butions by layer to the time mean difference in the CLR due to 
moisture differences. Between 50 and 90% of the major CLR 
differences arc due to pw differences between 200 and 500 hl\a 
and another 10-30% comes from the layer between 500 and 
700 hPa. 

We next extend the results of the previous analysis to con- 
sider the case where, in addition to the total CLR measure- 
ments, we also have available CLR measurements for the win- 


dow region. Figure 9 shows the sensitivity parameters f? 
computed separately for the window region (RadWn, le't-hand 
side panels) and for the total CLR minus RadWn (right-hand 
side panels). The latter consists of the rotation and vibration 
and rotation (RVR) bands or nonwindow bands; see disjussion 
of Figure 2V A comparison of Figures 9a and 9d she hat lha 
T tJ sensitivity in the window region tends to peak in middle 
latitudes, while in the nomvindow or RVR, the sensitivity 
shows a general increase with latitude. As such, the \ reatest 
sensitivity to 7 (J over the United Slates, Europe, and parts of 
China occurs in the window region. A comparison of figures 
9b and 9o with Figure 5a shows that the pw h sensi ivity is 
dominated by the RVR bands or nonwindow bands. Over the 
tropical oceans the sensitivity to the low-level moiMurc is 
greatest in the window region (compare Figures 9c and 9f). 

The large sensitivity over the eastern North Pacific and Atlan- 
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Figure 8. Relative contribution (%) of layer fipw to 5CLR. Contour intervals are for values opposite to Hie 
shaded values. 


lie seen in the total CLR calculation (Figure 6a) is due to the 
nonwindow bands or RVR bands (Figure 9f). 

By considering the linearization (2) for both the RadWn and 
the RVR bands f we obtain two equations relating the CLR 
differences in each band to changes in the geophysical param- 
eters. If we further limit the approximation of the CLR change 
to be primarily the result of just two geophysical parameters, 
we can estimate or “retrieve” these parameters. The previous 
results suggest that over land the CLR differences can be 
approximated by 

RadJFn ~ nr 7#(Rlt ^ (7) 

5RVR = ofr # <RVH)57' g + (^) 

while over ocean, the approximation is 

5 RadWn = ct rH ., m w n) ?>pWi + tvun.iuw ( 9 ) 

5RVR = <Xp»,iRVR)8p w i + il r~kmvn)fy nv h- (1^) 

Note that in the above relationships the sensitivity parameters 
(the ct M in Figure 9) arc obtained from the 6 hourly data 


according to (5). By inverting (7)-(8) and (9)-(I0), we obtain 
the retrieved estimates of the differences in the geophysical 
parameters assuming the and SRadWn and 5R V R arc 
known. For example, the retrieved monthly mean differences 
in the geophysical parameters ( T g , pw hf and/;^) arc shown 
m Figure 10 (left panels). The actual differences in T r , pw h , Flo 
an i\pw ( from the full radiative transfer equation shown earlier 
arc presented again in the right-hand panels to facilitate the 
comparison with the estimated values. The global mean bias in 
the retrieved 5T f/ is -O.O^C. The global mean root-mcan- 
Squnre error is I.77°C, which corresponds to about l n /< of the 
standard deviation of the true values. The global mean bias in 
the retrieval of 5/jiv,, is -0.1 mm. The rms error is 0.26 mm, 
which corresponds to about 8% of the standard deviation of 
the true values. Similarly, the global mean bias in the retrieval 
of fymy is 0.4 mm. The rms error is 2.15 mm, which corre- 
sponds to about 20% of the standard deviation of the true 
values. 

The above results show that (7) through (10) give a close 
approximation to the mean differences in the ground temper- 
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Figure 9. Coefficient of linear fit from Figures 9a lo 9c are based on .SRadVVn and from Figures W to 9f 
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mm" 1 ). 


ature and moisture and suggest that these relationships can be 
used to infer errors in the geophysical parameters. In the 
following section we examine the differences between the syn- 
thetic GHOS2 CLR and the CLR from a preliminary CERES/ 
TRMM data set. 

5. A Comparison With CERES/TRMM CLR 
and RadWn 

In the previous section we compared synthetic dear-sky 
fluxes (CLR and RadWn) computed from two 4-DDA data 
sets to develop some simple but quantitative measures of the 
sensitivity of clcar-sky fluxes to changes in ground temperature 
and moisture. Insight into the sensitivity of the clcar-sky fluxes 
to changes in these geophysical quantities was obtained by 


computing separately the various band contributions lo the 
clcar-sky fluxes. Reasonable approximations lo the sensitivity 
of clcar-sky fluxes to changes in the geophysical parameters 
were obtained with simple linear relationships using ground 
temperature and moisture in two broad layers of the atmo- 
sphere (the upper and lower troposphere) to predict th ; clcar- 
sky flux changes in the RadWn and RVR bands (or nonwindow 
bands): the single predictor relationships (5) and (6) as well as 
the two parameter relationships (7) through (10). 

We found that compared with ECMWF, GEOS2 Ins sub- 
stantially colder ground temperature during January 19' '8, with 
differences ranging from 3°C to more than 7°C over h gh lat- 
itudes. Over these regions the sensitivity of CLR t< is 
between 1.2 and L8 Win" 2 per PC, so differences in ll c CLR 
range from about 5 lo 15 Wm 2 in these regions, while parts of 
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the Antarctic show differences exceeding 20 Win” 2 . We also 
showed (Figures 10c and I Of) that compared with ECMWF, 
GEOS2 is welter in the upper troposphere during January 
1998 with differences as large as 2.5 mm over the tropical 
oceans. In the lower troposphere, GEOS2 is drier over most of 
the tropical oceans, with differences as large as 10 mm. 

The above results give us confidence in the methodology, 
though they only provide us with difference (or “error") fields 
relative to tire ECMWF products. We will now use some pre- 
liminary CERES/TRMM total CLR and RadWn data for Jan- 
uary 1998 (in place of the ECMWF analysis data) to evaluate 
the quality of the GEOS2 ground temperature over land and 
moisture over the oceans. 

The results for the moisture are shown in Figure 11. We will 
validate our results against the SSM/I total precipitable water 
(TPW). For purposes of illustration and to demonstrate the 
reliability of our algorithm, we first apply our methodology to 


the synthetic 5CLR and 5RadWn data. Here we assui ic that 
the ECMWF data arc the ground truth. The results Figure 
lie) show that indeed, we have adjusted the GEOS2TPW field 
to be close to that of ECMWF. 

We next demonstrate our algorithm by applying it to the 
satellite data (CLR and RadWn) and assume that the jalcllitc 
data are the ground truth. Figures 1 1 cl— Ilf show the results of 
applying our algorithm using the January 1998 CERES obser- 
vations of 5CL.R and 5RadWn. Figure I Id indicates liial our 
moisture between 200 and 500 hPa is excessive over the trop- 
ical regions. Figure lie indicates that at low levels GF.OS2 is 
too dry in subtropical regions around ±20°-30° in both hemi- 
spheres. The corrected GEOS2 TPW (based on Figures 1 td 
and He) is shown in Figure Ilf. Comparing with Figure I la, 
we see that the corrections have removed the major bias in the 
GEOS2 TPW fields. We note some amount of ovcrcori cclion, 
especially in the Northern Hemisphere subtropics. 
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Figure II, Diagnose of q: synthetic studies or TPW differences (mm): (a) GEOS2 minus SSM/I, (b) 
ECMWF minus SSM/I, and (c) GEOS2 corrected TPW based on ECMWF data minus SSM/I, arid satellite- 
based TPW and layer (d) retrieved &pw h based on ACLK and SRiulU'ii from CERES, (e) the same as Figure 
lid except for dpWf, and (f) GEOS2-corrcc!cd TPW based on Figures lid and 1 lc. Contour intervals arc for 
negative values, which arc opposite to the shaded values. 


Finally, we diagnose 7 f/ for January 199S. Since we do not 
have global ground temperature observations, we will again 
fi z compare with the ECMWF analyses. Figure 12a is the differ- 

ence map of the retrieved 5T f; minus actual (Figure 10a minus 
Figure lOd) based on the ECMWF data. As noted previously, 
the global mean bins is -0.04°C, and standard deviation is 
1.25°C. This again gives us confidence that our algorithm can 
correct the ground temperature to within those error limits. 
Figure 12b is the mean difference in the CLR between GEOS2 
and the CERES observations over the land areas between 4()°S 
and 40°N (the fraction of the globe covered by the CERES 
instrument). Differences arc negative over the Sudan, Saudi 
Arabia, and India. Positive values occur over southeastern 
Australia, southern Africa, the Andes, Central America, north- 
ern South America, and the Tibetan highlands. Figure 12c 


shows the implied GEOS2 T lf errors using (lie algorithm de- 
scribed above. The errors are quite large in some areas. For 
example, over the Sudan GEOS2 is more than 10°C too cold, 
while over southeast Australia, GEOS2 is more than I2°C too 
warm. Determining the errors in these estimates of the GEOS2 
ground temperature bias will require, among other thhgs, re- 
liable estimates of the bias in the CERES measurements. 

6. Discussion mul Conclusions 
We have introduced a simple but quantitative method for 
relating errors in model-based estimates of clcar-sky lo lgwavc 
(luxes (CLR and RadWn) lo errors in geophysical para neters. 
The primary motivation for this work is the underlying ? isump- 
lion that by linking the radiation errors to errors in t ic geo- 
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Figure 12. Diagnose of T tf \ synthetic studies of &T, f (°C) 
based on ECMWF: (a) GE0S2-rctricved 5 T tJ minus 5T, f hc- 
iwccn ECMWF and GEOS2. Satellite diagnose or AT,,': (b) 
8CLR between CERES and GEOS2, and (c) retrieved 57\ ; . 
Contour intervals are for negative values, which arc opposite 
to the shaded values. 


physical parameters, we can provide greater insight to the 
model developer on potential model errors. We show for a test 
ease (January 1998) that the method can potentially be used to 
obtain quantitative estimates of errors in ground temperature 
( <) T , f ) and moisture (fi/nv) from satellite observations. In par- 
ticular, our analysis of synthetic total and window region dear- 
sky ilux differences (computed from two different assimilated 
data sets) shows that a simple linear regression employing &T, t 
and broad layer Spw provides a good approximation to the full 
radiative transfer calculations, typically explaining more than 
90% of the 6 hourly variance in the flux differences. These 
simple regression relations can be inverted to "retrieve” the 
errors in the geophysical parameters. Uncertainties (normal- 
ized by standard deviation) in the monthly mean retrieved 
parameters range from 7% for 5 T y to about 20% for the lower 
tropospheric moisture 

Our initial application of the methodology employed an 
early CERES/TRMM data set (total and window region dear- 
sky longwave fluxes) to assess the quality of the GEOS2 data. 
The results showed that over the tropical and subtropical 
oceans, GEOS2 is, in genera!, too wet in the upper tropo- 


sphere (mean bias of 0.99 mm) and loo dry in the lower 
troposphere (mean bias of —4.7 nun). We note tha: these 
errors, as well as a cold bias in the ground temperntur;, have 
largely been corrected in the current version of GEO: 12 with 
the introduction of a land surface model, a moist turbulence 
scheme and the assimilation of SSM/I tola! prccipitabk water. 

The methodology described in this paper was developed as a 
result of our efforts to validate monthly mean Helds from the 
GEOS2 global data assimilation system, but the methodology 
should also prove useful for validating the climatological fields 
of global atmospheric models. The accuracy of the methodol- 
ogy depends on the accuracy of the radiation code, surface 
emissivity, the ozone profile, as well as the accuracy of the 
satellite estimates of total and window region dear-sky long- 
wave fluxes. While the initial results arc promising, further 
work is required to fully assess the sensitivity to error?, in the 
input parameters and to account for the mismatch between 
satellite observations and grid-scale model fields. 
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